Investigation of the dd — > 3 Hen7r° reaction 
with WASA-at-COSY 
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An exclusive measurement of the dd — > 3 Hen7r° reaction was carried at a beam momentum of 
Pd — 1.2 GeV/c using the WASA-at-COSY facility. For the first time data on the total cross section 
as well as differential distributions were obtained. The data are described with a phenomenological 
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approach based on a combination of a quasi-free model and a partial wave expansion for three-body 
reaction. The total cross section is found to be otot = (2.89±0.01 s tat. ±0.06 S! /s. ±0.29 nO rm) fJ,b. The 
contribution of the quasi-free processes (with the neutron being target or beam spectator) accounts 
for 38% of the total cross section and dominates the differential distributions in specific regions of 
the phase space. The remaining part of the cross section can be described within a partial wave 
decomposition indicating the significance of p-wave contributions in the final state. 

PACS numbers: 13.75.-n,21.45.-v,25.10.+s,25.45.-z 

Keywords: deuteron-deuteron interaction, pion production, few nucleon system 



INTRODUCTION 

At the fundamental level of the Standard Model, 
isospin violation is due to quark mass differences as well 
as electromagnetic effects [l|-|3j- Therefore, the observa- 
tion of isospin violation in principle allows one to study 
quark mass effects in hadronic processes. However, in 
general isospin violating observables are largely domi- 
nated by the pion mass differences, which are enhanced 
due to the small pion mass. An exception are charge sym- 
metry breaking (CSB) observables (charge symmetry is 
the invariance of a system under rotation by 180° around 
the second axis in isospin space that interchanges up and 
down quarks), since charge symmetry transforms a ir + 
into a 7r~ and therefore the pion mass difference does 
not contribute. Ref. [4| calls the investigation of CSB 
effects one of the most challenging subjects in hadron 
physics. On the basis of theoretical approaches with a 
direct connection to QCD, like lattice QCD and effective 
field theory, it is therefore possible to study quark mass 
effects on the hadronic level, since the effects of virtual 
photons are under control — for a detailed discussion on 
this subject see Ref. [5j. 

The first observation of the dd — > 4 He7r° reaction was 
reported for beam energies very close to the reaction 
threshold At the same time information on CSB in 
np — > dir° manifesting as a forward-backward asymme- 
try became available [7j. These data triggered advanced 
theoretical calculations within effective field theory, pro- 
viding the opportunity to investigate the influence of the 
quark masses in nuclear physics [8|, 0| . This is done us- 
ing Chiral Perturbation Theory (ChPT) which has been 
extended to pion production reactions First steps 

towards a theoretical understanding of the dd — > 4 Hc7r° 
reaction have been taken [ll, 12 1. Soft photon exchange 
in the initial state could significantly enhance the cross 
sections for dd — »■ 4 He7r° [13|, however, it was demon- 
strated in Ref. [3] that a simultaneous analysis of CSB 
in the two-nucleon sector and in dd — > 4 He7r° strongly 
constraints the calculations of the latter. 

The main problem in the calculation of dd — > 4 Hc7r° 
is to get theoretical control over the isospin symmetric 
part of the initial state interactions, for here high ac- 
curacy wave functions are needed for dd — > AN in low 
partial waves at relatively high energies. To get access 
to these a measurement is necessary for other, isospin 



conserving, dd induced pion production reaction chan- 
nels at a similar excess energy, such that the final state 
is constrained to small angular momenta. Then the in- 
coming system shares at least some of the partial waves 
in the initial state with the reaction dd — > 4 He7r°, while 
the transition operator is calculable with sufficient accu- 
racy using ChPT. Such a reaction is dd — > 3 Hen7r° and a 
corresponding measurement is presented here. 



EXPERIMENT 

The experiment was carried out at the Institute for Nu- 
clear Physics of Forschungszentrum Jiilich in Germany 
using the Cooler Synchrotron COSY [15| together with 
the WASA detection system. For the measurement of 
dd — > 3 Hen7r° at an excess energy of Q « 40 MeV a 
deuteron beam with a momentum of 1.2 GeV/c was scat- 
tered on frozen deuteron pellets provided by an internal 
pellet target. The reaction products 3 He and ir° were de- 
tected by the Forward Detector and the Central Detector 
of the WASA facility respectively, while the neutron re- 
mained undetected. The Forward Detector consists of 
several layers of plastic scintillators for particle identifi- 
cation and energy reconstruction and an array of straw 
tubes for precise tracking. The polar angular range be- 
tween 3° and 18° fully covers the angular range of the 
outgoing 3 He with the exception of very small angles. 
At this beam momentum the 3 He ejectiles have kinetic 
energies in the range of 65 - 214 MeV and, thus, are al- 
ready stopped in the first detector layers: in addition 
to the straw tube tracker only the two 3 mm thick lay- 
ers of the Forward Window Counter and the first 5 mm 
thick layer of the Forward Trigger Hodoscope were used. 
The two photons from the 7r° decay were detected by the 
Scintillator Electromagnetic Calorimeter as part of the 
Central Detector. The photons were distinguished from 
charged particles using the Plastic Scintillator Barrel lo- 
cated inside the calorimeter. The experiment trigger was 
based on a coincidence between a high energy deposit 
in both layers of the Forward Window Counter together 
with a veto condition on the first layer of the Forward 
Range Hodoscope to select helium ejectiles and a low en- 
ergy neutral cluster (E > 20 MeV) in the calorimeter to 
tag the decay of the pion. Further information on the 
WASA-at-COSY facility can be found in Ref. HI]. 
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FIG. 1. Left: Energy loss in the Forward Window Counter 
versus energy loss in the first layer of the Forward Trigger 
Hodoscope. The obtained energy pattern shows clear sepa- 
ration between different particles types. The graphical cut 
indicated in black represents the region used to select 3 He 
candidates. Right: The two photon invariant mass distribu- 
tion corresponding to the it — > 77 decay. The red dashed 
line indicates the tv° mass. 



DATA ANALYSIS 

Apart from the charge symmetry breaking reaction 
dd — > 4 He7r° with a four orders of magnitude smaller 
cross section, dd — >■ 3 Hen7r° is the only process with a 
charge 2 particle and a neutral pion in final state. Thus, 
the identification of a forward going helium and two neu- 
tral tracks forming a pion already provides a clean sig- 
nature for this reaction. Helium isotopes are identified 
by means of AE — AE plots using the energy deposit in 
the Forward Window Counter and the first layer of the 
Foward Trigger Hodoscope (Fig. [I] left). The condition 
that the 3 He has to pass at least the first two scintillator 
layers introduces an additional acceptance cut of Eki n > 
125 MeV. This rejects most of the 3 He going backward 
in the cm. system. However, having two identical par- 
ticles in the initial state and, thus, a symmetric angular 
distribution with respect to 9cm = 90° the full angu- 
lar range can be recovered by a symmetrization of the 
detected events. The energy deposits are also used to 
reconstruct the 3 He kinetic energy by matching the en- 
ergy loss pattern to Monte-Carlo simulations. The 3 He 
four-momentum is completed by the direction informa- 
tion from the straw tube tracker. In addition to the 3 He 
two neutral clusters in the central detector corresponding 
to the two photons from the 7r° decay were requested. As 
event pile-up and low energy satellites of genuine photon 
clusters can cause larger photon multiplicities the most 
probable true two-photon combination was identified by 
selecting the pair with the 3 He — ir° missing mass being 
closest to the neutron mass. As result a nearly back- 
ground free pion peak was obtained (Fig.[l] right). In a fi- 
nal step the data were refined by applying a kinematic fit 
using the hypothesis dd — > 3 Hen7r°. Still remaining back- 



ground and badly reconstructed events were rejected by 
a cut on the cumulated probability distribution at 10%. 
At the end of the analysis chain about 3.4 x 10 6 fully 
reconstructed and background free dd — > 3 Hen7r° events 
are available. It should be noted that the Dalitz plot is 
fully covered except for a small region for large 7r° — n in- 
variant masses due to the acceptance hole for #3h c < 3°. 
Although based on this data set any possible differen- 
tial distribution can be generated — e.g. for a selective 
comparison with future microscopic theoretical calcula- 
tions — a suitable set of observables for further analysis 
and presentation had to be selected. For any unpolarised 
measurement with three particles in final state four inde- 
pendent variables fully describe the reaction kinematics. 
For the present analysis the choice for these independent 
variables is based on the Jacobi momenta q and p with 
q being the ir° momentum in the overall cm. system 
and p the momentum in the rest frame of the 3 He — n 
subsystem. The following variables were constructed ac- 
cordingly: cos6 qi cosOp (the polar angles of q and p, re- 
spectively), M3 Hc „ and ip (the angle between the projec- 
tions of q and p onto the xy- plane) . As discussed earlier 
all plots show data after a symmetrization in the global 
cm. system. 

The absolute normalisation was done relative to the 
dd — > 3 Hen reaction. Corresponding data were taken in 
parallel during the first part of the run using a separate 
trigger. Due to the correlation between kinetic energy 
and scattering angle for the binary reaction, quasi mono- 
energetic particles form a distinct and clean peak in the 
AE — AE plots. For the selected events the 3 He missing 
mass distribution reveals a background free peak at the 
mass of the neutron (Fig. [3J left). In order to determine 
the integrated luminosity the data presented in Ref. (l7| 
were used. The authors measured the reaction dd — > 3 Hj> 
for beam momenta between 1.09 GeV/c - 1.78 GeV/c and 
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FIG. 2. Measurement of the dd — » 3 Hen reaction. Left: 
3 He missing mass distribution, the vertical red dashed line 
indicates the neutron mass. Right: Measured angular distri- 
bution in comparison with a Monte Carlo (MC) simulation. 
The MC sample was generated based on a parametrized cross 
section. The cross section is given in arbitrary units (a.u.). 
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dd — > 3 Hen for beam momenta in the range of 1.65 GeV/c 
- 2.5 GeV/c. Moreover, they showed that the differential 
cross sections for both channels at 1.65 GeV/c are iden- 
tical within the presented errors. Based on these results 
we used the measured cross sections for dd — > 3 Hp to 
calculate the cross sections for dd — > 3 Hen at 1.2 GeV/c. 
For this the angular distributions for the beam momenta 
of 1.109 GeV/c, 1.387 GeV/c and 1.493 GeV/c were 
parametrized. Then, for selected polar angles the de- 
pendence of the differential cross section on the beam 
momentum was fitted and interpolated to the beam mo- 
mentum of 1.2 GeV/c. The resulting distribution was 
used as an input for the simulation of dd — > 3 Hen. Fig- 
ure [2] (right) shows the match of the angular distribution 
of 3 He in data and the Monte-Carlo filtered event gener- 
ator. The extracted integrated luminosity is determined 



t0 be L \nt 



(877 ± 2 s tat. ± 62 sys . ± 62„ orm .) nb~ 



where the superscript 1 refers to the first part of the 
run, see below. The systematic uncertainty reflects dif- 
ferent parametrizations of the reference cross sections. 
In addition, the uncertainty of 7% in the absolute nor- 
malization of the reference data is also included. The 
result for the total cross section given below is based 
only on this first part of the run. The second part was 
optimized for high luminosities and also served as a pi- 
lot run for a measurement of dd — > 4 He7r°. It provided 
data to extract high statistics differential distributions 
for dd — > 3 Hcn7T°. These have been absolutely normal- 
ized relative to the first part of the run using the rates 
of the dd — > 3 Hen7r° reaction. The integrated luminos- 
ity obtained for the second part of the run amounts to 
L 2 mt = (4909 ± 13 stai . ± 35Q sys . ± 350„ orm .) nfe- 1 . 

The uncertainty on the integrated luminosity (in total 
10% if all contributions are added quadratically) is the 
dominant source for the systematic error on the abso- 
lute normalisation. Another major source is associated 
with the cut on the cumulated probability distribution 
of the kinematic fit. In order to quantify the influence of 
this cut, the analysis was repeated for different regions 
in the probability distribution. For the total cross sec- 
tion the maximum deviation from the average value was 
taken as error. Changes in the shape of the differential 
distributions were extracted similarly, however excluding 
the variation in the absolute scale. For all other analysis 
conditions according to the criteria discussed in Ref. [18j 
no significant systematic effect was observed. 



PHENOMENOLOGICAL MODELS 



based on dp 3 He7r° and a partial- wave expansion for 
the 3-body reaction. While the latter is limited to s- 
and p-waves, the large relative momentum between the 
spectator nucleon and the rest system in the quasi-free 
model corresponds to higher partial waves motivating the 
incoherent sum and the neglection of interference terms. 



Quasi-free reaction model 

High momentum transfer reactions involving a 
deuteron can proceed via interaction with a single nu- 
cleon of the deuteron with the second nucleon being re- 
garded as a spectator. Naturally, this mechanism is most 
significant in regions of the phase space where the mo- 
mentum of one nucleon in final state matches the typical 
Fermi momenta in the deuteron. In the present experi- 
ment two deuterons are involved and, thus, the reaction 
may proceed with a projectile or target neutron spec- 
tator. For the parametrization of the quasi-free sub re- 
action dp — > 3 He7r°, the empirical angular distributions 
and the energy dependent cross section in the energy 
regime from threshold up to an excess energy of 10 MeV 
[20j and for excess energies of 40, 60 and 80 MeV 21] 
have been used. They have been convoluted with the 
momentum distribution of the proton in the deuteron 
using an analytical form of the deuteron wave function 
based on the Paris potential [l9(. As a result one gets 
absolutely normalized differential cross sections for the 
quasi-free contribution to dd — > 3 Hen7r°, which can be 
directly compared to the measured data. Figure [3] (left) 
shows the momentum distribution of the neutron for data 
and the quasi-free model filtered by Monte-Carlo. As ex- 
pected the quasi-free process dominates the distribution 
for small momenta. The lower boundary of the spectrum 
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Presently, no theoretical calculation exists for a mi- 
croscopic description of the investigated reaction. How- 
ever, in order to have a sufficiently precise acceptance 
correction a model which reproduces the experimental 
data reasonably well is required. The ansatz used here 
is the incoherent sum of a quasi-free reaction mechanism 



FIG. 3. Comparison of data (black line) and the quasi-free 
model filtered by Monte-Carlo (blue dashed line). Left: Mo- 
mentum distribution of the neutron. Right: Angular distri- 
bution of the pion in the 3 He — n° subsystem for neutron 
momenta smaller than 90 MeV/c (indicated by the vertical 
red dashed line in the left plot). Data are not corrected for 
the acceptance. 
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is caused by kinematic effects. At a beam momentum 
of 1.2 GeV/c the reaction dp — > 3 He7r° with the tar- 
get proton at rest is below threshold and can only oc- 
cur for pfermi > 48 MeV/c. The vanishing acceptance at 
#3 Hc < 3° further increases the minimum Fermi momen- 
tum. The right panel of Fig. [3] shows the angular distri- 
bution of the pion in the 3 He — ir° subsystem for neutron 
momenta below 90 MeV/c, i.e. in the region where the 
quasi-free process should dominate the distribution. 

Partial wave decomposition 

For the remaining part of the data which cannot be de- 
scribed with a quasi-free process a 3-body model based 
on a partial wave decomposition has been developed. 

I 



2tt dM3n cn d cos 9 p dcos9 q dip L 

i A 2 q 2 (1 + 3 cos 26 q ) + |Atp 2 (1 + 3 cos 26 p ) 



The relative angular momenta were defined according to 
the coordinates introduced earlier: one in the global ir° 
- ( 3 Hen) system and one in the 3 Hen subsystem (de- 
noted by I and L, respectively). For the partial wave de- 
composition the angular momenta have been limited to 
I + L < 1, i.e. to at most one p-wave in the system. For 
the momentum dependence the standard approximation 
\M\ 2 cx q 2l p 2L was used. Taking into account all possible 
spin configurations this results in 18 possible amplitudes. 
After combining the amplitudes with the same signature 
in final state, four possible contributions can be identi- 
fied: s-wave in both systems (sS), one p-wave in either 
system (sP and pS) and a sP — pS interference term. 
They can be described by seven real coefficients (four 
complex amplitudes minus one overall phase). With this 
the four-fold differential cross section can be written as: 



A iq 2 + A 3 p 2 + 

(1) 

f A 5 pq cos 9 p cos 9 q + A 6 pq sin 9 p sin 6 q cos ip } 



ith 



C = 



32(27r) 5 sp*(2s a + l)(2s 6 + l) 



(2) 



where s a and s& denote the spin of beam and target and 
s and p* a the total energy squared and the beam momen- 
tum, respectively, in the cm. system. The coefficients 
Ai describe the strength of the individual contributions 
mentioned above: A$ corresponds to I = L = (sS), A-y 
and A 2 to I = 1 and L = (pS), A 3 and Ai to I = 
and L = 1 (sP) and A$ and Aq to the interference term. 
Integration of Eq. [1] results in a set of equations for the 
description of the single differential cross sections: 



da 



dM3 Kcn 

da 
27rdcos 9 n 



16n 2 C P q [A + A iq 2 + A 3 p 2 ] (3a) 



da 



27r<icos 9 V 



da 

dp 



AttC 
AttC 
8irC 



B + -A 2 (1 + 3 cos 26» q )I pS 
B + jAi (1 + 3cos26» p )/ s p 



B + — A 6 I pS+s p cos ip 
16 



(3b) 
(3c) 
(3d) 



with the new coefficient 



B = A I sS + AJ p s + A 3 I sP . 



(4) 



The constants 7 S 5, I p s, I s p and I p s+ S p are the results of 
integration over M3 Hen : 



I P s — 



hp — 



tpS+sP 



(M 3Ho +Af„)2 

(M 3Ho +Af„)2 
(x/i-Af,) 2 

(A/ 3Ho +Af„) 2 
(v^-Af x ) 2 



pqdM 3lic 



pq 3 dA'h Hc 



p 3 qdA'h Rc 



p 2 q 2 dM 3lic 



(5a) 
(5b) 
(5c) 
(5d) 



'(A/ 3Ho +Af„) 2 

Equations [3] do not contain the coefficient A$ as the 
corresponding term vanishes with the integration over 
cosO q and cos9 p . In order to extract this coefficient Eq. [1] 
has to be multiplied by cos9 q cos9 p before integration. 
This results in the following formula to determine A 5 : 



da ' 8 . 

— = -irCA 5 I pS+sP 



(G) 



with a'(q,p) — a(q,p) ■ cos9 q cos9 p . 

It has to be noted that the coefficients Aq, Ai and ^3 
cannot be extracted unambiguously from the differential 
distribution da/dMa Hcn . In the non-relativistic limit q 2 
and p 2 are both linear in M3^ cn introducing a correlation 
of all three coefficients. For the measurement of dd —> 
3 Hen7r° at an excess energy of Q ss 40 MeV this is still 
a good approximation. Thus, only a value for B can be 
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FIG. 4. Acceptance corrected data (black points) presented 
as functions of M3jj e nj cos8 q , cos8 p and cp, respectively. The 
curves represent the fit to the model: full model (red solid), 
quasi-free contribution (blue dashed) and the partial wave de- 
composition (green long dashed). The hatched areas indicate 
the systematic uncertainties of the shape of differential distri- 
butions. Uncertainties on the absolute normalization are not 
included. 



extracted from the data. Any values for Aq, A\ and A3 
fulfilling Eq. |4]and the fit to da/dMs-Ren W1U to the 
same model description. 



RESULTS 

In a first step a sum of Monte-Carlo filtered distribu- 
tions for each contribution from the partial wave decom- 
position (coefficients Aq to Aq) and from the quasi- free 
model (coefficient A-f) was fitted to the uncorrected, sin- 
gle differential spectra. The result served as input for the 
Monte-Carlo simulation finally used to determine the ac- 
ceptance correction. 

The final distributions after acceptance correction are 
presented in Fig. 2) Contributions from the quasi-free 
model, the partial wave decomposition and the full model 
are shown in blue, green and red, respectively. These 
spectra were refitted using the analytical formulas given 
in the previous section. The result is consistent with the 
previous fit. Although the partial wave expansion was 
limited to at most one p-wave in the final state it pro- 
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5. Distribution of da' /dip as used in Eq. [6] to extract 



vides a reasonable overall description of the data: both 
angular distributions show a significant contribution of 
p- waves of similar size, the pS-sP interference term is vi- 
sualized by the non- isotropic distribution of da /dip. The 
quasi-free contribution is about 1.11 \xb and, thus, is in 
agreement within the normalisation error of about 10% 
given in Ref. [2(| with the prediction of the quasi-free 
model of 1.192 The result of the fit using Eq. |6]and 
the quasi-free model to the data is presented in Fig. [5] 
The values for extracted coefficients from the global fit 
are summarized in Table HI 

One should emphasize that the meaning of the fit pa- 
rameters is limited to the context of the discussed model. 
Any addition of higher partial waves, for example, would 
also change the extracted amplitudes of the lower par- 
tial waves. Thus, systematic errors are only provided for 
data. For the extracted fit parameters only statistical 
errors are given. 



Contribution 


Fit results 


sS,pS,Sp 


B [lib] 


wave 


1.840 ± 0.003 


sS 


A [fib/GeV] 


wave 


(0.41 ± 0.08) xlO 4 


pS 


Ai [^6/ GeV& ] 


A 2 [fib/GeV 11 ] 


wave 


(8.2 ± 7.0) xlO 4 


(18.3 ± 0.3)xl0 4 


sP 


A 3 [^b/GeV°] 


A 4 [^/ GeV& ] 


wave 


(1.08 ± 0.02) xlO 4 


(18.04 ± 0.07) xlO 4 


pS+Sp 


A 5 ht&/GeV & ] 


A 6 [fib/GeV] 


interference 


(-45.4 ± 0.3) xlO 4 


(-15.0 ± 0.2) xlO 4 


quasi free 


CT q f ■ Ar [fib] 




1.108 ± 0.003 



TABLE I. Collection of the extracted fit parameters. The 
amplitudes are given in units of (47r) 2 C. The parameters Ao, 
Ai and ^3 are correlated and could not be extracted unam- 
biguously, the given numbers represent one possible solution 
(see text). 



So far, a possible momentum dependence of the tran- 
sition amplitudes, for example due to initial or final state 
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He-n excess energy / MeV 



differential distributions are only one possible represen- 
tation of the results. One goal of the measurement was to 
provide data for studying dd initial state interaction for 
small angular momenta, which is one missing informa- 
tion in the microscopic description of the charge symme- 
try breaking reaction dd — > 4 He7r° within the framework 
of Chiral Perturbation Theory. Once the important ob- 
servables have been identified the corresponding experi- 
mental distributions can be provided. 
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FIG. 6. Coefficient A4 representing the p-wave contribution 
in the 3 He — n system as a function of the excess energy. The 
error bars along the x-axis represent the width of the intervals 
in M 3Hen . 

interaction, was neglected. Deviations from this assump- 
tion were studied by refitting the data for five intervals 
in M3fjen (corresponding to intervals in q and p). All co- 
efficients except one remained constant. Only A4 repre- 
senting a p-wave contribution in the 3 Hen system showed 
a significant momentum dependence (see Fig. [5]): A4 is 
larger for low excess energies in the 3 Hert system (corre- 
sponding to low relative momenta). One possible reason 
for this might be excited states with isospin I = 1 in 
the 3 Hen system at low excess energies as reported in 
Ref. (the production of an I = state would be 
charge symmetry breaking). 

Integrating over the differential distributions we obtain 
for the total cross section of the dd — > 3 Hen7r° reaction: 

a tot = (2.89 ± 0.0 W ± 0.06 sys . ± 0.29„ orm ) fib. (7) 
SUMMARY 

For the first time an exclusive measurement of the 
dd — > 3 Hen7T° reaction has been performed. A total cross 
section of o tot — 2.89 fib with an accuracy of about 11% 
has been extracted. Differential distributions have been 
compared to the incoherent sum of a quasi-free reaction 
model and a partial-wave expansion limited to at most 
one p-wave in the final state. The contribution of the 
quasi- free processes accounts for about 1.11 fib of the to- 
tal cross section matching the prediction of the quasi-free 
reaction model. The partial wave decomposition reveals 
the importance of p-wave contributions in the final state. 
The applied model shows a reasonable agreement for all 
differential distribution. Thus, based on this compari- 
son no indication for significant contributions of higher 
partial waves can be deduced. 

The whole data set amounts to about 3.4 x 10 6 fully re- 
constructed and background-free events. The presented 
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